Matter stability in modified teleparallel 

gravity 



A. Behboodi, S. Akhshabi, K. Nozari 

Department of Physics, Faculty of Basic Sciences, University of 
Mazandaran, P. O. Box 47416-9544 7, Babolsar, IRAN 



Abstract 

We study the matter stability in modified teleparallel gravity or f(T) theo- 
ries. We show that there is no Dolgov-Kawasaki instability in these types of 
modified teleparallel gravity theories. This gives the f(T) theories a great 
advantage over their f(R) counterparts because from the stability point of 
view there isn't any limit on the form of functions that can be chosen. 
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1. Introduction 

It was Einstein who soon after formulating his theory of general relativity, 
first introduced the idea of teleparallel gravity [0]. In this new theory, a set 
of four tetrad (or vierbein) fields form the orthogonal bases for the tangent 
space at each point of spacetime and torsion instead of curvature describes 
gravitational interactions. Tetrads are the dynamical variables and play a 
similar role to the metric tensor field in general relativity. Teleparallel gravity 
also uses the curvature-free Weitzenbock connection instead of Levi-Civita 
connenction of general relativity to define covariant derivatives j^j. 

After its first introduction, further important developments were made by 
several pioneering works in teleparallel gravity and it has been shown that 
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teleparallel Lagrangian density only differs with Ricci scalar by a total diver- 
gence 0, 0] • This shows that general relativity and teleparallel gravity are 
dynamically equivalent theories where the difference arises only in boundary 
terms. However there are some fundamental conceptual differences between 
teleparallel theory and general relativity. According to general relativity, 
gravity curves the spacetime and shapes the geometry. In teleparallel theory 
however torsion does not shape the geometry but instead acts as a force. 
This means that there are no geodesies equations in teleparallel gravity but 
there are force equations much like the Lorentz force in electrodynamics. 

Recently with the discovery of accelerated cosmic expansion modifying 
gravity beyond general relativity has generated much interest. One way to 
modify gravity is to replace the GR Lagrangian density, R, with a general 
function of Ricci scalar. This approach leads to the so called f(R) theories 
of gravity 0]. Similarly one can try to modify gravity in the context of 
teleparallel formalism and replace the teleparallel Lagrangian density, T with 
a general function of T which leads to the generalized teleparallel gravity or 
f(T) theories Q. The resulting field equations in f(T) theories are second 
order equations and are much simpler than the fourth order equations that 
appear in metric formalism of f(R) gravities. 

It has been shown that f(T) theories can explain the present time cosmic 
acceleration without resorting to some exotic dark energy @, 0, 10]. However 



one should remain cautious when selecting the form of function f(T). It is a 
well established fact in our every day experience that weak-field gravitational 
bodies like the sun or the earth do not experience violent instabilities result- 
ing in dramatic changes in their gravitational fields. So any theory which 
results in such instabilities should be clearly ruled out. It has been shown 



that some prototypes of f(R) theory suffer from these instabilities [ll| and 
a general condition for the stability of such theories has been derived 12 . 
Similarly one should consider stability of the theory in the weak-field limit 
of f(T) gravity. In this paper we show that matter is generally stable in the 
context of modified teleparallel gravity. 

2. Field equations 

In teleparallel gravity we need to define four orthogonal vector fields 
named tetrad which form the basis of spacetime. The manifold and the 
Minkowski metrics are related as 

9^ = Vije'ei (1) 
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where the Greek indices run from to 3 in coordinate basis of the mani- 
fold, the Latin indices run the same in tangent space of the manifold and 
rjij = diag(+l, —1, —1, —1). The connection in teleparallel theory, the 
Weitzenbock connection, is defined as 

IV = ef^e; (2) 

which gives the spacetime a nonzero torsion but zero curvature in contrast 
with general relativity. By this definition the torsion tensor and its permu- 
tations are \^ 

T% u = e?(d„el - d v ei) (3) 
K"% = -l(T^ p -T^ p -T p n (4) 

sr = \{K»\ + - %T°* a ). (5) 

Where is called the superpotential. In correspondence with Ricci scalar 
we define a torsion scalar as 

T = S/»T^ U (6) 

so the gravitational action is 

where |e| is the determinant of the vierbein which is equal to \J—g. Varia- 
tion of the above action with respect to the vierbeins will give the teleparallel 
field equations 

e-%(ee?S p n - e'T^S^ + \e\T = AirGe^ (8) 

Now similar to modifying the action of general relativity which R is replaced 
by a general function f(R), one can replace the teleparallel action T by a 
function f(T). Doing this, the resulting modified field equations are 

e-%(eS/ v )f(T)-e^ x S p ^f(T) 

+ S { ^d,{T)f"{T) + \e»f{T) = 47rGef 6; (9) 

where Q p is the energy momentum tensor of matter. In what follows we set 
AnG = 1. 
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3. Matter Stability 

The main motivation for modifying gravity in both teleparallel and gen- 
eral relativity is the explanation of present time accelerated expansion of the 
universe. If one considers a flat, homogeneous Friedmann- Robertson- Walker 
universe, then the tetrads are 

= diag(l, a(t), a(t), a(t)j (10) 

and the torsion scalar will be 

T = -6^ = -6H 2 (11) 
a 2 

From the field equation (9) one can derive the modified Friedmann equation 
as 

12H 2 f'(T) + f(T) = 4p (12) 

To achieve the present time acceleration, any added term to the torsion scalar 
should be dominant at late times but negligible at early times. In ref. [§[ 
the form f(T) = T — e/(— T) n has been proposed. This gives the correct 
cosmological dynamics at late times without resorting to dark energy. 

Now we turn to the problem of matter stability. Following the above 
discussion we promote the torsion scalar, T to a general function in the form 

f(T)=T + e<p(T) (13) 

where the parameter e should be small to agree with recent observational 
constraints. To study the matter stability of a model of modified teleparallel 
gravity in the form (13), we begin by taking the trace of field equation (9) 

+ Tf{T) + S^d,(T)f"(T) + f(T) = (14) 
Substituting (13), gives 

+T(1 + e<ft) 

+ S v ^d,(TW +(T + ep{T)) = 6 (15) 
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Note that Eqs. (14) and (15) correspond to the trace of the equation of 
motion since the only kind of perturbations we take into account are the 
ones of conformal factor i.e. scalar modes. Now we apply this equation to 
the gravitational field of a weak-field object like the sun or the earth. For 
such gravitational bodies, the torsion scalar in linear perturbation can be 



approximated by [11 



- T = 6 + 2V M T/ I/ + T x (16) 

where T\ is the linear perturbation and V M is the covariant derivative with 
the Levi-Civita connection. This equation followed from the fact that the 
torsion scalar T and the Ricci scalar R only differ by a total divergence, 
R = —T — 2V ' pT^ 11 . The minus sign in (16) comes from the fact that 
the torsion scalar is negative for a homogeneous and isotropic weak field 
gravitational body. 

The metric is also approximately can be taken as the Minkowski metric 
plus some small perturbations 

9HU = TJfMu + h^y (17) 

where we assume perturbations to be homogeneous and isotropic. Eq. (17) 
means that the vierbeins can also locally be written in the form 

e: = 6j + ej (18) 

where ej is a small perturbation in relation to the trivial tetrad. We can 
describe the deviation from the flat spacetime by [4] 

ej = ae( 1)v + a 2 e{ 2)u + ... (19) 

where a is a dimensionless parameter which labels the order of perturbations. 
Inserting this expansion in (17), the corresponding expansion of the metric 
is 

9p.v = Vim + <x(e(i)vw + + ... (20) 

and we have efa u = S-e^ and e^v = V^(\) v - 

Here we consider only the first order or linear perturbations so from now on 
we drop the subscript (1) from the equations. 

By perturbing the torsion in the form of equation (16) we'll have the following 
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equation in linear perturbation theory for the nearly flat region inside a weak 
field celestial body (see Appendix for proof) 



Tl0 + Ve^K 



e 



2Ve 



(21) 



where A = is a positive constant. Inserting (16) , (18) and (21) in 
(15) and keeping only the terms linear in perturbations, we get 



A Aecp' . „ r- 
H %= + Aecp'WQ 



2ve 2ve 



-A(l+€(f') 



+ 




- Ae(tp' - cp") 







+ etp' + 2 



40V© 



40V0- 

-A(l+etp') \d t (e)+Ve df ~e v v \ +A(2+e V ') 

A 







L2V0 



d 2 e v 



-A{1 + e<p')Ved t (e»)5l + -V0ey/0 

0eo9 eo9 

2 Y 4 ^ 



(22) 



where A = and a dot denotes differentiation with respect to time. Note 
that the perturbation equation in modified teleparallel gravity, equation (22) 
is a first order differential equation in contrast to the second order equations 
that appear in f(R) theories 12] . The right hand side of (22) is a source term 
involving the matter content and also deviation from the flat background as 
in (17) and (18). Equation (22) can be rewritten in a concise form as 



where we have defined 



m 



mTi + nTi = IT 



(23) 



2V0 2V© 



n = -A(l + e<p') 



n ee — Ae(<p' — if") 








40V0 



40V0 
A(l + e<p' 



+ e<p' + 2 
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+A{2 + eip') 



e 



Qd t (e v 



e&p" + e - ^ee^' 



-etp 



(24) 



-2ve 

Let us make a comparison between values of the terms in m. For a typical 
gravitational body the energy- momentum scalar, G, is proportional to the 
mass density of the body and is positive [11] 



~ (10 3 sec)- 2 ( 



Pr. 



g cm 



(25) 



where p m is the mass density of the body. For example we have p m = 
5.52g/cm 3 for the earth and p m = 1.41 g/cm 3 for the sun. The value of 
e is fixed in such a way that it gives the correct cosmological dynamics at 
late times, so it should be extremely small. For example, a common class of 
functions that are popular in f(T) literature is 



f(T) = T — 



2(n+l) 



CO' 



(26) 



where n is some real number and the p parameter will be fixed to a value that 
the model can reproduce the late time accelerated expansion of the universe 
[§, [h|. For this model we have 



10 18 sec 



(27) 



From this it is obvious that the first term in m is much larger than the 
other two terms and we can safely neglect the second and third terms. Doing 
this, equation (22) becomes 



,_^„ )( e ) + ^e + ^ 



2VQ 

A 



n 



(28) 



Let's consider the time evolution of perturbations. From the form of 
differential equation (28) it is obvious that first order perturbations, T\ will 
grow with time if the coefficient of 7\ in (28) is negative and decreases with 
time if the coefficient is positive. Growing of perturbations with time will 
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mean that the torsion will rise very quickly and leads to strong instability 
while a decreasing perturbations will mean that the gravitational field will 
bounce back to its equilibrium state and so the body is stable. The coefficient 
of Ti in (28) is dominated by the last term due to extremely small value 
of e. Note that A and 6 are positive so from this discussion it is obvious that 
the coefficient of Ti will always remain positive and as a result the matter in 
these types of theories is always stable. 

Now we turn our attention to the case of a radiation fluid. For this type 
of matter the trace of the energy- momentum tensor, is vanishing. From 
(16) we have — T = V^T/" + T\ = T{ . Inserting this in the trace of field 
equation (9) yields 

f[ + p Tf /2 - q T{ 1/2 = (29) 

where by definition 

P = 




Solving equation (29) for the time evolution of Ti gives 

T[ (t) = ^tanh (^ty/pq+ ^Vpo)" (31) 

which of course is always stable because the perturbations will become con- 
stant after some time. Here C is an integration constant. The limiting value 
is given by q/p = Aeip/ (2 + eip) which is extremely small because of the value 
of e. Figure (1) shows the qualitative behaviour of Ti as given by equation 
(31). 



4. Conclusion 

From a geometric point of view, modifying gravity seems a necessary task 
in order to explain recent positively accelerated expansion of the universe. 
Any such modified theory, whether it is in the context of general relativity 
or in teleparallel gravity, may be expected to show some strong deviation 
from the standard gravity at very high energies and in strong-field regimes. 
This is because we still do not have a proper theory of quantum gravity to 
describe the behavior of gravitational interactions at those energies. On the 
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Figure 1: Qualitative behavior of the first order torsion scalar perturbation versus 
time for radiation matter with vanishing energy- momentum scalar. T{ will reach 
a constant value given by q/p and so the matter in this scenario is always stable. 

other hand any strong deviation from the standard gravity at low energies 
and weak-field regimes immediately disqualify the theory because it will con- 
tradicts well established weak-field experiments. One of these experiments 
is the stability of weak-field celestial bodies or any other weak gravity ob- 
jects. In this paper we've investigated the stability of such objects in the 
context of modified teleparallel gravity. The analysis shows that there is no 
Dolgov-Kawasaki matter instability in these type of theories. In contrast, in 
the corresponding f(R) theories a certain stability condition should be met. 
This gives a great advantage to f(T) theories over their f(R) counterparts 
because from matter stability viewpoint, there is no limit on the form of 
functions that can be chosen to replace the torsion scalar in the action of 
f(T) theories. We note that we have extended our analysis to the second 
order of perturbations and we have observed that the matter is still stable in 
this scenario. 
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Appendix A. 

Here we present the proof of equation (23) for an almost flat region inside 
a weak field gravitational body. For such an object the tetrad and metric are 
given by equations (18) and (20) respectively. Considering only the first order 
perturbations and dropping the subscript we have the following equations for 
the torsion and superpotential tensors 

T' lw = d lt ef-d v e£ (A.l) 

and 

S/f = d p e» - d»e p - - d a e ap ) 

+ S P(d^;-d a e^) (A.2) 

the tensor is not necessarily symmetric but it has been shown that the 
anti-symmetric part of it has no physical significance in the field equations so 
we assume it to be symmetric here [4]. Furthermore for an almost flat region 
inside a star, we can safely assume that both the background and the first 
order correction are homogeneous and isotropic. In that case the torsion and 
its perturbation does not depend on spatial coordinates and we have <9 M — > d t . 
Also for a homogeneous and isotropic perturbation the first order correction 
of the tetrad has the form 

= diag(l,b,b,b) (A.3) 

and b only depends on time. Substituting this in (Al) and (A2), we can find 
the torsion scalar as 

T = S ppv T p ^ = -6b 2 (A.4) 

Up to the first order in perturbations, the second term in the R.H.S of Eq. 
(16) will be V = 36. On the other hand the only non zero components 
of the superpotential tensor are all the same (up to a sign) and proportional 
to y/T, in particular we have 

d,Sr = h (A.5) 
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so we will have the relation 

d,Sr = ^jft^) ( A -6) 
substituting from (16), equation (23) is obtained. 
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